Percolation of particles on recursive lattices: 
(II) the effect of size and shape disparities 

Andrea Corsi*, P. D. Gujrati''" 
The Department of Physics and The Department of Polymer Science, 
The University of Akron, Akron, Ohio, 44325, U.S.A. 
(Dated: February 6, 2008) 

Abstract 

The preparation of many composites requires the intermixing of several macromolecular fluids 
along with the addition of rigid filler particles. These fillers are usually polydisperse and there is 
an extensive experimental evidence that their size and shape profoundly affects the properties of 
the resulting material. In particular, it is generally found that the percolation threshold decreases 
as the size disparity between the different particles present in a system increases and that the 
threshold decreases with the aspect ratio of the particles. Here, a recursive approach that we 
have recently introduced is applied to the study of the percolation of particles of different sizes 
and shapes, without the presence of a polymer matrix, on a lattice in various phases including 
metastable states. In our approach, the original lattice is replaced by a recursive structure on which 
calculations are done exactly and interactions as well as size and shape disparities are easily taken 
into account. In the previous paper of this series, we have introduced the recursive approach and 
shown how correlations among particles of the same size can affect percolation. Before considering 
the complete system made of particles of various sizes and shapes embedded in a polymer matrix, 
in the third paper of the series, we describe here the properties of systems made of particles without 
any matrix. The approach appears to be extremely successful since it is able to capture most of 
the important features observed in experiments. 
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I. INTRODUCTION 



The fabrication of typical polymeric systems often requires the intermixing of several 
macromolecular fluids along with the addition of rigid filler particles. The most classical 
example of inclusion of rigid particles in a polymer matrix is represented by the addition 
of carbon black into a natural rubber matrix in order to improve its strength and process- 
ability. The original driving force for the incorporation of fillers into polymers was certainly 
economic since common polymers produced in high volume typically cost about ten times 
more than mineral fillers. In most cases, though, the cost argument cannot be sustained. 
The incorporation of fillers into thermoplastic polymers through compounding, for example, 
is a costly process that almost always cancels out the economic advantage of introducing 
the filler. 

The use of powdered materials of sizes ranging from nanometers to microns is not limited 
to polymer composites but is common in many branches of materials science. Fillers are 
carefully selected to enhance the performance of a matrix. Fillers have been used in order 
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properties of the host material. The elastic modulus of a filled resin, as an example, results 
from a complex interplay of the properties of the individual components. There is enough 
experimental evidence to conclude that the properties of the composite are affected by a large 
number of parameters: the size, shape and distribution of the reinforcing particles as well as 
the interactions between the particles and the polymeric matrix. Even without considering 
the interactions of the particles with the polymer matrix, it is clear from the experiments 
that the percolation threshold decreases as the size disparity between the different particles 
present in a system increases and that the threshold decreases with the aspect ratio of 
the particles present in the system J, Q, ^ 19, Is^. In the case of particles that are not 
spherical in shape, another important parameter that affects the mechanical properties of 
the composite is represented by the degree of orientation of the particles with respect to the 
direction of an applied external stress state. 

Although most of the physical systems that are usually described through percolative 
models are polydisperse, most studies have described the percolation of equal sized com- 
ponents. Even though lattice models are simpler, the effect of polydispersity has been 
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considered only in continuum percolation, mainly in the case of a system of particles repre- 
sented by discs and spheres with a distribution of radii. In contrast to some initial claims 
that the average fractional volume covered by the percolating species at the threshold would 
be constant 26|, results that have been recently published have shown that the volume 
fraction occupied by the percolating species at the threshold increases as the ratio between 
the size of the constituent species becomes larger and larger, more significantly so when the 
ratio itself is very large [^l • 

Many other models have been developed in order to describe the percolation process. 
Some of them involve modifications of the classical theory of rubber elasticity to account for 



additional cross-links produced by the fillers 



Models have also been developed based 



on a double network description j29| and on the replica formalism 



Monte-Carlo lattice simulations of filled 
on the rotational isomeric-state model ^ 
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yrners are also very common and often based 



3^,|33i. The simulations are applied to both 



polymer chains attached at one end to the filler surface and to free chains. All these sim- 
ulations assume that the filler particles occupy the sites of a regular cubic lattice and the 
polymer chains occupy the space between the particles. These simulations were successful 
in predicting the anisotropic reinforcing properties of prolate particles 



II. OBJECTIVE OF THE PRESENT RESEARCH 



As explained before [3J|, a composite that contains nanoparticles is characterized by 
nanoscopic inhomogeneities that cannot be described satisfactorily in a continuum model 
since the averaging processes performed during the calculations are done on length scales 
much larger than the size of the nanoparticles. A lattice model is more promising and can be 
used to capture this nanoscopic inhomogeneities. Previous lattice model calculations of per- 
colation generally did not take into account the possible difference between the size of filler 
particles (with very few exceptions like 3^) and usually neglected interactions. Such calcu- 
lations have additional very important limitations represented by the use of random mixing 
approximation, the incompressibility of the model, and the necessity for the monomers and 
voids to have the same size. 

In our approach, we replace the original regular lattice by a recursive lattice (RL), which 
is built up from its smaller parts in a recursive fashion. The two infinitely large recursive 
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structures that we have extensively used are the Bethe lattice and the Husimi lattice. The 
choice of the recursive lattice to be used is dictated by the model being investigated. 

Our theory is given by the solution on a recursive lattice. These tree structures allow 
us to capture only weak correlations and are consequently not suitable, for example, for 
carrying out the calculation of critical exponents. However, other kinds of recursive lattices 
that do not have a tree structure can be used to obtain non-classical exponents, as we have 



done in 



36|. However, as our interest here is not in critical exponents, we will restrict 



ourselves to the theory obtained on tree-like recursive lattices. This theory has been applied 
by Gujrati and coworkers to study and describe a wide range of polymer systems and it has 
provided extensive insight in their phase separation, critical points, loop formation in tree 



polymer gel, theta states, compressibility effects, immiscib 
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Here we study the possible percolation of filler particles of different sizes and shapes in the 
system. Although size disparity alone without any interactions has been argued to induce 
phase separation in this kind of system, recent rigorous calculations [s^ have proven that 
no phase separation in an athermal fully packed state of hard particle mixtures on a lattice 
is possible merely due to size disparity. 

Our goal is to study the dependence of percolation on the size disparitv and nature of 
the interactions between different filler particles. In the following paper we describe 
the effects of the presence of a polymer matrix on the percolation of these particles. 

Even though conventionally, when a solution is described, the name solvent is attributed 
to the abundant species, in the following we will always refer to the smallest component 
present in the system as the solvent even when that species is not the abundant one. This 
does not mean that different species must be chemically different, it is just a convention to 
make the description easier to understand. 

We study a system of monomeric and square particles and compare its behavior with a 
system of monomeric and star particles on a square Husimi lattice to determine the effects of 
the change in shape of the larger particles on the properties of systems made of particles of 
different sizes. In this and the following publications, we will simply call the square Husimi 
lattice the Husimi lattice for brevity. 
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III. MONOMERIC AND SQUARE PARTICLES 



Model 



The first system we investigate is the one formed by monomeric particles and square 
particles on a Husimi lattice, see Figures ^ and |21 Figure ^ shows a portion of a Husimi 
lattice, including the labeling of its sites. The Husimi lattice is built in a recursive fashion 
starting from a central plaquette. The different size of the plaquettes that belong to different 
generations is just a matter of convention: one might consider all plaquettes to have the 
same size but then a lattice would become impossible to draw because of the crowding that 
one has as one moves away from the origin of the lattice. Figure |2] shows one element of a 
Husimi lattice, that is a square plaquette. As shown in the figure, the square particles that 
we are considering are particles in the shape of squares occupying one elemental square of 



the Husimi lattice and, therefore, its four lattice sites. As explained before j3J|, we define 
"above" and "below" relative to the origin of the lattice (above meaning moving away from 
the origin) and we distinguish the sites inside a plaquette of the Husimi lattice as a base site 
(the one closest to the origin of the lattice), two intermediate sites and a peak site, opposite 
to the base one. 

In order to be able to describe this system, we introduce three possible states that describe 
the configuration of the mth level site which is the base site of the mth level square. This 
site can be in the S state if it is occupied by a solvent molecule, in the A state if a corner of 
a square particle occupies the site and the square is above it (so that all the (m + l)th level 
sites are occupied by the remaining corners of the square) or in the B state if a corner of a 
square particle occupies the site and the square lies below it, see Figure El 

As we have seen in the first part 3^, we can associate a Boltzmann weight w with every 
nearest-neighbor contact between particles of different species present in the system. The 
weight w is determined by the excess interaction energy e as follows: 

w = exp {—/3e) . (1) 

We will set e =1 to set the temperature scale, as was also done in the first paper. The 
activity for the monomeric species particle, is r], related to its chemical potential /x through 

r] = exp (/9/i) = w'''. (2) 
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The total partition function for this system of sites can be written as: 

Z = 5^ti;%^= (3) 

where Nc is the number of contacts between particles of two different species and A'^^ is the 
number of solvent particles present in the lattice. From the partition function of the system 
we can obtain the thermodynamic potential or free energy per site F by means of 

Z = e-^^/^^ (4) 

In the following, instead of using the conventional expression of the free energy (per site), 
we will use the dimensionless free energy 

00 = -pF. (5) 

The change in sign implies that whenever we will look for stable phases, we will have to look 
for the maximum of the free energy and not for its minimum. Taking into account, as it 
usually happens, that the sum in the expression for Z over all sets of states is determined by 
the equilibrium state corresponding to the average densities per site 0c and 0s, where (pc is 
the density of contacts between solvent molecules and square particles and 0s is the density 
of solvent molecules, and given the definition of the entropy of the system, we can write for 
the entropy per site 

S =UJ - /?/i0s + /30c. (6) 

All thermodynamic densities in this work will be defined per site unless noted otherwise. The 
procedure that we use to obtain the above two densities has been outlined in the previous 
paper 3^: it is necessary to consider the contribution to the total partition function given 
by configurations in which the solvent is at the origin (to obtain the value of 0s) or by those 
in which the particle at the origin is in "contact" with one or more particles of a different 
species (to obtain the value of 0c). 

B. Recursion Relations 

When the mth level site is in the S state, there are eight possible configurations for the 
three sites at the (m + l)th level as shown in Figure El 
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Following js^, we can write the recursion relation for Zm{S), the partial partition function 
for the mth branch of the Husimi lattice given that the mth level site is occupied by a solvent 
molecule, as a function of the Zm+i- 

(S) = r][Zm.+i (S) (S) (S) 
+ 3w^Zm+i (S) Zm.+i (S) (A) 
+ {2w' + m;^)Z„,+i (S) Zm+i (A) Zm+i (A) 
+ w^Z„,+, (A) Zm+i (A) Zm+i ( A)] . (7) 

If the the mth level site is in the B state, we can write a similar expression for Z^ (B) 
since the possible configurations of the three sites at the (m + l)th level are the same as in 
the previous case: 

Zm (B) = w'^Zm+1 (S) Zm.+l (S) Zm+1 (S) 

+ i2w^ + w')Zm+l (S) Zm+1 (S) Zm+1 (A) 
+ 3w'^Zm+l (S) Zm+1 (A) Zm+1 (A) 

+ Zm+1 (A) (A) Zm+1 (A) , (8) 

where the Boltzmann weight w appears every time a solvent molecule is a nearest neighbor 
of a square particle. The possible configurations of the (m + l)th level sites when the base 
site is in the B state are shown in Figure |3] Finally, if the the mth level is in the A state 
we can write: 

Zm (A) = Zm+1 (B) Zm+1 (B) Zm+1 (B) . (9) 

The single configuration that the system can assume in this particular case is shown in 
Figure 

We introduce the following ratios: 

XmiS)^^,Xm{A)^^^ (10) 



and 



where we have introduced 



Xm (B) ^ = 1 - (S) - Xm (A) , (11) 



Bm = Zm (S) + Zm (A) + Zm (B) . (12) 
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The ratios satisfy the following sum rule at every generation 

XmiA)+XmiB)+XmiS) = l. (13) 

We express Bm in terms of the partial partition functions of the (m + l)th level for com- 
pleteness: 

Bm = Vl^m+l (S) Zm+l (S) Zm+1 (S) + 3w'^Zm+l (S) Zm+1 (S) Zm+1 (A) 

+ i2w^ + W^)Zm+l (S) Zm+1 (A) Zm+1 (A) + W^Zm+1 (A) Zm+1 (A) Zm+1 (A)] 

(S) Zm+1 (S) Z^+i (S) + {2w' + W )Zm+l (S) (S) (A) 
+ 3w^Zm+i (S) Z^+l (A) Zm+1 (A) + (A) Zm+1 (A) (A) 

+ Zm+1 (B) (B) (B) , (14) 

which is easy to derive. 

The choice of the normalization factor Bm used to obtain these ratios from the partial 
partition functions is critical. We always try to keep these ratios finite so that the denomi- 
nator cannot go to zero. In order to enforce this, we always use combinations such as that 
above: our lattice is always completely filled with either squares or solvent molecules present 
at any temperature. Thus, at any temperature, the terms Zm (S), Zm (A) or Zm (B) cannot 
all be zero and all the ratios are always finite. At the fixed point of the recursion relations 
we have: 

Xm (S) = Xm+l (S) = S, Xm (A) = Xm+1 (A) = O, Xm (B) = Xm+1 (B) = b. (15) 

In general, we call such a fixed point solution, in which the ratios do not change from 
level to level, a 1-cycle solution. The reason for this nomenclature will become clear below. 
Using the fixed point solution, we can express 

Bm = Bl^.Q, (16) 

where Bm+i = Zm+i {S) + Zm.+i (A) + Zm+i (B) according to equation[T2l and where we have 
introduced a polynomial Q, which at the fix-point takes the value 

Q = {ri + w^)s'^ + {3r]w^ + 2w^ + w^)s^a + {3w^ + 2riw'^ + riw^)sa'^ + (1 + w^)a^ + b\ (17) 
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Using this polynomial and the recursion relations written above, we can write at the fix-point 
the following system of equations to be solved: 



sQ = r][s'^ + 3w^s^a + (2w^ + w^)sa^ + w'^a% (18a) 
aQ = 6^ (18b) 
bQ = wh^ + {2w^ + w'^)s^a + Sw^sa"^ + a^. (18c) 



C. Phase Diagram and Percolation 



1. Densities 



In order to obtain the densities of the two species, we must consider the total partition 
function of our system. The total partition function of the system at the (m = 0)th level can 



be written considering all the possible configurations of the (m = 0)th level site j5^. The 
entire lattice is obtained by joining two (m = 0)th level branches Co; see for definition. 
In order to obtain the total partition function, we must consider all the configurations that 
the system can assume in the two branches that meet at the origin of the lattice. 

There are only three possible configurations that the system can assume at the origin, as 
shown in Figure El The site at the origin can be occupied by either a solvent molecule or the 
corner of a square particle. In the case of a square particle, we have to take into account the 
fact that the square can be either above or below the origin. The total partition function of 
the system can then be written as 

^^^Zo(S)Zo(S)^2Zo(A)Zo(B). (19) 
V 

The first term represents the contribution of the configuration with the solvent molecule at 
the origin, while the second one represents the contribution of the configurations with the 



square at the origin. As explained before [3j], the factor rj in the denominator of the first 
term is necessary in order not to over-count the solvent activity. The factor two in the last 
term is related to the possibility of having the square on either side of the origin. 

The form of this contribution is easily understood if we realize that every configuration 
that appears as an A state looking at it from one side of the origin, appears to be in the B 
state if considered from the opposite side. 
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The solvent density is then defined as the ratio between that part of the total partition 
function containing configurations in which the origin is occupied by a solvent molecule and 
the total partition function of the system. 

In this particular case, we can express the density of the solvent particles on the lattice 

as: 

Z, (S) Zo (S) h 

4>s = ^ • (20) 

Similarly, the density of sites covered by the square particles is defined as 

2Zo (A) Zq (B) 

K = ^ . (21) 

This is, by the definition, the mass density of squares. It is obvious that 0sq + 0s = 1- Since 
each square has four corners (and occupies four sites), the number density of squares is 



0sq/4. (22) 



2. Free Energy per Site 

The total partition function can be used to obtain the thermodynamics of the system. 
It is clear that Zq is the total partition function of the system obtained by joining two 
branches Co together at the origin. In the thermodynamic limit, Zq becomes unbounded 
and care must be exercised. For this purpose, we need to consider the free energy per site. 
This we accomplish as follows. Let us imagine taking away from the lattice the two squares 
that meet at the origin. This leaves behind six different branches Ci. We connect these 
branches to form three smaller but identical cacti, the partition function of each of which is 
denoted by Zi; the latter can be written in a form that is identical to that of equation 
except that the index of each partial partition function is replaced by 1. 

The difference between the free energy of the complete lattice and that of the three 
reduced lattices is just the free energy corresponding to a pair of squares (which contain 4 



sites of the lattice) so that, following Gujrati 



52| | , we can write the adimensional free energy 



per site without the conventional minus sign as 

1 1 / Zq \ 

u = -(cj (complete lattice) — 3u;(reduced lattice)) = ^ i^J ' ^^"^^ 

It is possible to write 

Zo = B^oQ2,Zi = BlQ2, (24) 
10 



where Bm is defined above and Q2 is, at the fixed point, the polynomial 



Q2 = — + 2a6. 
V 



(25) 



Since Bq = BfQ, see (fTBj). the free energy per site can be written as 



LU = 




(26) 



The calculations are done in the grand canonical ensemble and are carried out at constant 
chemical potentials. Hence, the above free energy represents iSPgqVo where /5 is the inverse 
temperature. „„ represents the volume of the uuit eel. of the lattice aud P., is the os,not,e 
pressure |53( across a membrane permeable to the square particles. The osmotic pressure 
Ps across a membrane permeable to the solvent particles is given by (3PsVq = (3PsqVQ — Inr]. 
If the solvents represent voids, then Pg represents the conventional pressure of the lattice 
system js^. The approach that we use is the following: for every value of w (T) we solve 
the system of recursive equations that we have derived. Then we use the fixed point values 
of a, b and s to calculate the value of the free energy at that w. 

In order to determine which phase is the stable one at some temperature, we must find 
the free energy of all the possible phases of the system as a function of w. 



Let us first consider the case when /i is positive. The ground state of the lattice at 
zero temperature is represented by a pure solvent phase in this case: all the lattice sites 
are occupied by monomeric species. As we start increasing the temperature, the density of 
squares increases. In this case, the solvent is percolating at every temperature. Thus, its 
percolation is of no interest. 

For positive values of fi, only one solution of the set of recursive equations exists and 
consequently only one phase is present and, therefore, no phase transition is observed. From 
what was said earlier, this phase is easily seen to be a mixture of solvent particles and squares, 
which contains only solvent particles at absolute zero. As no phase transition occurs, we 
need to only investigate is the percolation of squares to which we now turn. 

To study the percolation of the square particles, we introduce as before the probability 
-Rm(A) < 1 that a site in an A state at the mth generation is connected to a finite cluster of 



3. /i > 
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squares at higher generations as well as a probability i?m(B) < 1 that a site in a B state at 
the mth generation is connected to a finite cluster of squares at higher generations. Then 
Zjn{oi)Rm{(y) denotes the contribution to the partial partition function of the branch Cm due 
to all those configurations in which the site at the mth generation is in the state a and is 
connected to a finite cluster of squares at higher generations. If we divide by 
Zm{ot) we obtain a recursion relation for the Rm{ot)- At the fixed-point solution, each one 
of these Rm{o) approaches its fixed-point value given by the solution of an equation of the 
form 

R{a)=p[{Rmi (27) 

where {R{(3)} represents the set formed by the complete set of states (3. By considering the 
configurations shown in Figures ID and El and following an approach that is identical to that 
used to obtain the partial partition functions, we can write: 

R{k) = i?(B)3, (28) 
w^s^ + {2w^R{A) + w^)s^a + 3w^R{Afsa^ + RjAfa^ 

This system of equations does not depend explicitly on the chemical potential fi but it 
depends on it implicitly. In fact, for every value of T, and hence of w;, a different chemical 
potential corresponds to different values of s, a and b and consequently a different value for 
the R{a)'s. 

Since R{a) represents the probability of the (occupied) site at the origin of the lattice 
being connected to a finite cluster of squares on one half of the lattice and since two half 
lattices are joined together at the origin to form the complete lattice, we define 

Psq = 1 - R{A)R{B) (30) 

as the percolation probability for the squares. This represents the probability of the occupied 
site at the origin being connected to an infinite cluster of squares on both sides of the origin. 
The behavior of psq as a function of 0sq for different (positive) values of fi, is shown in Figure 
ini As expected, the percolation threshold increases as the chemical potential increases. 
The athermal case, when the two species do not have any interactions, corresponds to 
fi = +00, e — > such that fie > is kept constant. 

As before, we can study the percolation process as a function of the temperature. Figure 
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[7| shows the percolation probabihty as a function of the temperature for the same chemical 
potential values considered in Figure IHl 

Figure |H1 shows how the percolation threshold depends on the value of the chemical 
potential while Figure IHl shows the dependence of the critical temperature. As expected, an 
increase in the chemical potential for the smaller particles makes it more and more difficult 
for the bigger particles to percolate. This explains the increase in the critical temperature 
and the percolation threshold as /i increases. As before, the limiting value for the percolation 
threshold as fi goes to infinity is the athermal value while the critical temperature grows 
unbounded. 

4. ^l<o 

If negative values of fi are considered, the problem becomes more interesting and slightly 
more complicated. A negative value for fi means that, according to thermodynamics, the 
system will be stable at low temperature if no solvent is present in the lattice and all the sites 
are covered by square particles. The very nature of the square particles, in particular the fact 
that they occupy more than one single site on the lattice, makes the solution of the problem 
more challenging. What can be expected at low temperature, or at very low concentrations 
of solvent molecules is an ordered structure of squares. The difference between the nature of 
this phase and that of a pure solvent phase must be appreciated. The order present in the 
case of a pure solvent phase is due to the regular nature of the lattice (the periodic nature 
in the case of a regular lattice), and has nothing to do with the nature of the interactions 
between the particles that are present in the system. Instead, when a phase that is very rich 
in square particles is present, the order of the structure is associated with the nature of the 
particles and is not only due to the lattice itself. 

In the present case, we expect to have an ordered phase at low temperature that eventually 
disorders at high temperature, when large fractions of both solvent and square particles are 
present in the system. Thus, we expect to have a phase transition between an ordered 
phase, stable at low temperatures, and a disordered phase, stable at higher temperatures. 
The empty circles in Figure ^1 represent the free energy of the system as obtained by using 
the above method of 1-cycle calculation in the case of fi = —1. This cannot be the free 
energy corresponding to the stable phase at every temperature, since it becomes negative. 
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The adimensional free energy uj is the logarithm of the partition function. If we are properly 
describing the ground state of our system, the partition function at T = reduces to 2, since 
the only configuration that contributes is the lattice fully packed by square particles. On 
a square lattice, there are only two such configurations corresponding to the checkerboard 
packing of the lattice. Thus, the partition function at any positive temperatures contains 
this term besides other positive terms. Hence, lo for a stable phase cannot be negative. 
This means that the phase described by the 1-cycle solution is not describing the stable 
phase at low temperatures because of the phase transition. It is, therefore, necessary to 
introduce a different description of our configurations on the lattice in order to be able to 
capture the phase that is stable at low temperatures. This is a problem that we have already 



encountered in previous studies 
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43l ] and it requires a different fixed point solution than 



the conventional 1-cycle solution used for dealing with disordered phases. 

The presence of large square particles introduces a correlation between sites that are not 
nearest neighbors. In order to describe these correlations, let us first consider a system 
made of species that are all monomeric in size. If the base site of a particular mth level 
plaquette is occupied by a particle of one species, then the middle sites and the peak site 
of the same plaquette can be occupied in principle by a particle belonging to any species. 
In the presence of larger particles that occupy multiple sites of the lattice, this is not true 
anymore. A first example of this has been presented above when we were looking at the 
possible configurations of the system when the site of interest is in the A state, as shown 
in Figure El In this case the two middle sites and the peak site do not have any freedom 
and can only be in the B state. The recursion relation for Zm(A) in fact contains one single 
term. But the correlations introduced by the presence of these large particles affect even 
larger portions of the lattice. If we move one level beyond to the (m + l)th plaquette, we 
immediately understand that the presence of the square particle inside the mth plaquette 
affects the state of the sites in the (m + l)th plaquette as well. The middle sites and the peak 
site in the (m + l)th plaquette cannot be in the B state because that would be incompatible 
with the state of the sites one level below. From this discussion, it should be clear how 
the introduction of non-monomeric species induces correlations that would have not been 
present if we just had monomeric species. 

From the point of view of the recursion relations, it is necessary to look for a different 
kind of fixed point solution. When using our recursive approach, we always start from some 
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initial guesses for the ratios that describe the possible states of any site of the lattice. These 
initial guesses can be thought of as representing some particular boundary conditions at the 
surface of the system. Then, by using the appropriate recursion relations for the problem 
under investigation, we descend the lattice and we obtain the values of the ratios at the mth 
generation of the lattice as a function of the values at the (m + l)th generation. When the 
difference between the values of the ratios at two consecutive levels is lower than a preset 
tolerance, usually set equal to 10~^^, it is assumed that the fixed point has been reached and 
that the values obtained for the ratios describe the behavior in the bulk of the system. The 
values of the ratios at the fixed point can be used to obtain the free energy of the system 
and the entire thermodynamics as explained in the previous section. Different initial guesses 
are chosen in order to investigate if different fixed points, and consequently different phases, 
are present in the system for any temperature. This method, as explained above, is what 
we refer to as the one-cycle method. This explains the meaning of the label for the lower 
free energy curve in Figure ITnl 

This approach never captures the ordered phase of square particles like the one that is 
expected at low temperatures in the case of negative chemical potentials. To obtain the 
appropriate description for the ordered phase at low temperatures, we start with the state 
at absolute zero, where all the lattice sites are occupied by the square particles. If the base 
site (index m) of a plaquette is in the A state then the middle sites (index m + 1) and the 
peak site (index m + 1) of that plaquette are in the B state. The middle and peak sites 
(index m + 2) in the following plaquette are in the A state and so on. It is possible to notice 
a periodic structure in which consecutive generations on the lattice are alternatively in the 
A and B states. The one-cycle approach to the solution of the recursion relations will not 
succeed in describing this ordered phase. What we have is this: if at one generation m of 
the lattice a = 0, for example, then at the following generation m + 1, a = 1 followed by 
a = and so on so forth. In this case, the difference am — am+i never goes to zero, but the 
difference am — CLm+2 does. This provides us with the clue to describing the low temperature 
phase ordered phase where the system is mostly filled with square particles as follows. We 
need to use the two-cycle description for the ordered phase at low temperatures. In this 
case what is checked is the difference of the values of the ratios at two levels that are one 
generation apart on the lattice. Whenever — am+2 (as well as hm — &m+2 and Sm — Sm+i) 
but not am — Om+i (as well as hm — bm+i and Sm — Sm+i) is less than the desired tolerance. 
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we assume the fixed point of tfie solutions lias been reached. The system of equations for 
this kind of fix-point is the same as derived above and the equations are used to calculate 
the quantities of interest needed to describe the thermodynamics of the problem as well as 
the percolation probabilities for the system. 

It should be stressed that the two-cycle fix-point contains the one-cycle fix-point as a 
special case when the quantities at two successive generations become the same. Of course 
one might also obtain higher order cycles in which quantities such with n > 2 

go to zero as the fixed point is approached. For an example of such a situation, we refer the 
reader to the following section. In the case of the system under investigation we have tried 
to look at three-, four- and five-cycle solutions without getting any physically meaningful 
results. But we do have to point out that, at least in principle, one can not rule out the 
presence of higher order cycles without actually looking for them, even though they are hard 
to envision. One could expect a 4-cycle structure in the present case if there were attractive 
interactions {w > 1) such as specific interactions like hydrogen bonding between the solvent 
and the solute. But we do not consider these interactions. However, in general, the available 
cycle structures would emerge from the same set of recursive relations by starting from all 
possible initial guesses, and we have searched for different ones but we have been able to 
find just the two possible solutions mentioned in this section. 

By using this second method, it is possible to obtain a second fixed point solution. The 
corresponding free energy curve is the upper one represented by the circles Figure ^1 The 
solutions obtained with the two different schemes coincide at high temperatures, but are 
different at low temperatures. The temperature at which the ordered phase appears and 
becomes the stable solution is, in the case of /x = —1, Tqd — 15.7, indicated by the dotted 
line in the figure. We use the subscript OD to remind that this transition is a phase transition 
between an ordered and a disordered phase. It is possible to observe that the stable phase 
at every temperature has a positive free energy, as expected. Since u lacks the conventional 
minus sign, the stable phase is the one that has the maximum free energy. 

The entropy corresponding to the free energy shown in Figure El is shown in Figure ^2 

As expected, below Tqd the entropy of the ordered phase is lower than the entropy of 
the metastable continuation of the disordered phase. The entropy of the metastable phase, 
though, drops much more rapidly, and goes to zero at a finite temperature. This temperature 
is the Kauzmann temperature of the system, corresponding to the Kauzmann catastrophe 
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5J| . This result is very interesting since until recently the Kauzmann catastrophe had been 



observed theoretically only in polymeric systems |41|,|43|, but never in system made of small 
molecules. Recently Semerianov [55[ has shown the presence of a Kauzmann catastrophe in 
a dimer model. 

The results for the percolation probability as a function of the temperature are shown 
in Figure ^| It is obvious that the strength of the percolation process is very different 
for the two phases. The formation of a percolating cluster of solvent molecules is much 
easier in the disordered phase than in the ordered one. There is a wide temperature range, 
between ~ 8.42 and ~ 2.25 in the case fi = —1, in which there is at least one percolating 
cluster in the disordered phase but no percolation occurs in the ordered phase. The region 
below 2.25, limited by the dotted line in Figure IT^ does not represent physical states of the 
disordered phase because it is below the Kauzmann temperature and it corresponds to a 
negative entropy. 

The percolation probability is plotted as a function of the solvent density in Figure CHI 
This figure shows how the percolation threshold for the two phases is very different. In 
particular the solvent density must be at least ~ 0.34 in order to have percolation in the 
ordered phase while it is much lower in the metastable prosecution of the disordered phase. 
While lowering the temperature, if the phase transition is avoided the system is percolating 
all the way down to the Kauzmann temperature. 

Figures and El summarize the results obtained as the value of the solvent chemical 
potential is changed. As expected, an increase in the magnitude of the chemical potential 
makes it more and more difficult for the solvent particles to percolate because such an 
increase makes it more and more favorable for the lattice to be covered by square particles 
reducing as much as possible the number of solvent particles present in the system, at least 
at low temperature. 

If we consider the athermal limit in which there is no interaction between the two species 
of particles, we obtain the value (ps — 0.23 and (psq — 0.255. These results are very interesting 
for two reasons. First of all, we observe that the percolation threshold in this athermal case 
is larger for the square particles than for the solvent particles. This is consistent with most 
experimental findings that show how the percolation threshold increases with the size of 
the particles. Also, the percolation threshold is lower for both species with respect to the 
value obtained for same size species in other studies, (ps — 0.38. This behavior is also 
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consistent with most of the experimental findings that the presence of size disparity lowers 
the percolation threshold of all the species present in the system. 

IV. MONOMERIC AND 5-SITE STAR PARTICLES 
A. Model 

The next system that we have analyzed is the one formed by monomeric particles and 
star particles on a Husimi lattice, see Figure ITHl compare it with Figure |21 We only consider 
the smallest possible stars, each of which occupies five sites of the lattice, as compared to 
the four sites occupied by the square particles considered above, and always have their core 
(the center) on a corner of a square plaquette of the lattice. Consequently, the star particles 
always occupy sites belonging to two different and adjacent plaquettes of the lattice. The 
introduction of the star shaped particles introduces correlations among sites of the lattice 
that have a range that is longer than the one characteristic of the squares. This makes the 
problem a bit more challenging. In order to be able to describe this system we introduced 
four possible states that describe the configuration of the mth level site of the lattice. The 
base site in the mth level square can be in an A state if an end point of a star particle 
occupies the site and the core of the star is above it, a B state if an end point of the star 
occupies the site and the core of the star lies below it or a C state if the core of a star occupies 
the site, see Figure E| It is important to point out that the A and B states describe all 
configurations in which an end point is at the site of interest with a core above or below, 
respectively. So parts (a) and (b) of Figure El only show one of the two possible A and B 
configurations each. The star could be also on the other side of the plaquette as well. The 
site is in the S state if it is occupied by a solvent molecule. As before, we always refer to 
above and below with respect to the origin of the lattice: above means farther away from 
the origin of the lattice while below means closer to the origin. 

Following what has been done before in the case of square particles, we can associate 
a Boltzmann weight w with every contact between particles of different nature present in 
the system. Like before, the weight w is determined by the excess interaction energy e as 
w = exp {—/3e). The activity for the solvent, rj, is obtained as a function of its chemical 
potential, fi, through 77 = exp (Pfi) = as before; see 1^. 



18 



Similar to the case of monomeric and square particles, the total partition function for 
this system can be written as 0- The entropy per site is similarly given by (jH)). 

B. Recursion Relations 

When a site at the mth level is in the S state, there are nine possible configurations 
for the three sites at the (m + l)th level. These nine configurations are shown in Figure 
IT71 In the figure, the stars are always shown on one side but it is important to remember 
that, as explained above, they could be on either side of the plaquette. The configurations 
are very similar to those obtained in the previous section for the percolation of monomeric 
and square particles with the important difference that in the case of the stars it is possible 
to have the configuration in which the core of a star is at the peak site of the (m + l)th 
plaquette of the lattice with two arms of the stars being inside such plaquette (configuration 
(i) in the figure). In all configurations, the stars occupies sites of the lattice belonging to 
two different plaquettes. 

By using the same considerations as in the previous section, we can construct the recursion 
relation for Zm{S), the partial partition function for the mth branch of the Husimi lattice 
given that the mth level site is occupied by a solvent molecule as: 

Zm (S) = ?7[^.m+l (S) ^m+l (S) Z^+l (S) + 3w^Z^rn+l (S) Zm+1 (S) Z^+l (A) 

+ i2w^ + (S) Zm+i (A) (A) + w^Z^^.,, (A) Z^+i (A) Z^+i (A) 

+ (C) (B) Z^+i (B)], (31) 

where rj is the activity. 

If the site at the mth level is in the A state, we can write: 

Zm (A) = 2Zm+l (C) Zm+l (B) Zm+1 (A) + 2w'Zm+i (C) Zm+1 (B) Zm+1 (S) . (32) 

The configurations that the system can assume in this particular case are shown in Figure 

dl 

If the site at the mth level is in the B state, we can write an expression that is very 
similar to the one written above for the S state since the possible configurations of the 
three sites at the (m + l)th level are the same as in that case and only the weights have 
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to be properly changed. Then by properly taking into account all the nearest neighbors 
interactions between particles of different species, we finally have: 

(B) = w'Zm^^ (S) (S) (S) + i^w" + (S) Z^+i (S) Z^+i (A) 

+ 'iw'Zm^^ (S) (A) (A) + Z^+i (A) Zm+i (A) Z^+i (A) (33) 
(B)Z„+i(B), (34) 

where the Boltzmann weight w appears every time a solvent molecule is a nearest neighbor 
of a star particle. The possible configurations of the (m + l)th level sites are shown in Figure 
IT^ Finally, if the site at the mth level is in the C state, we have: 

Z^ (C) = Zm.+i (B) Zrn+l (B) Zm+1 (A) + W^Z^+i (B) Zm+1 (B) Zm+1 (S) . (35) 

The configurations that the system can assume in this particular case are shown in Figure 
M 

Following the approach introduced in the previous section, we introduce the following 
ratios: 

Xn,{S) = — — ,x„(A) = — ,x„(B) = — (36) 

-Om -Dm -Om 

(C) = = 1 - x„ (S) - (A) - x„ (B) . (37) 

-Dm 

where we have introduced 

= Z.^ (S) + Z.^ (A) + (B) + Zm (C) . (38) 
The ratios satisfy the following sum rule at every generation 

Xm (A) + Xm (B) + Xm (C) + Xm (S) = l] (39) 

compare with (fT^ . 

At the (1-cycle) fixed point of the recursion relations, we have: 



Xm (S) 


Xm+1 


(S) 


= s, 


Xm (A) 


Xm+1 


(A) 


= a, 


Xm (B) 


Xm+1 


(B) 


= b, 


Xm (C) 


Xm+1 


(C) 


= c, 



(40) 
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Then, at this fix-point, we can introduce a polynomial Q, as done in p6p . which is given by 



+ (1 + M;2)a^ + + ivw^ + l)cb^ + 2(a + m;2s)6c + (a + w^s)^^. (41) 

Using this polynomial and the recursion relations written above, we can write 

sQ = r][s^ + Sw^s^a + {2w^ + w^)sa'^ + w'^a^ + w^cb% (42) 

aQ = 2{a + w^s)bc, (43) 

bQ = w^s^ + {2w^ + w^)s'^a + 3w^sa^ + + cb'^, (44) 

cQ = {a + wh)b^ (45) 



C. Phase Diagram and Percolation 

1. Densities and Free Energy per Site 

In order to obtain the densities of the two species, we must consider the total partition 
function of our system. 

There are only three possible configurations that the system can assume at the origin, as 
shown in Figure ^1 The site at the origin can be occupied by either a solvent molecule, or 
the end point of a star or the core of a star. In the case of a star particle, we have to take 
into account the fact that the end point of the star can be either above or below the origin. 
The total partition function of the system can then be written as 

Zo = + 2Zo (A) Zo (B) + Zo (C) (C) . (46) 

The first term represents the contribution of the configuration with the solvent molecule 
at the origin, while the second one represents the contribution of the two possible config- 
urations with the end point at the origin and the last one represents the contribution of 
the configuration with the core at the origin. The factor two is related to the possibility of 
having the star on either side of the origin when the end point is at the origin. 

The solvent density is then defined as the ratio between the partition function describing 
configurations in which the origin is occupied by a solvent molecule and the total partition 
function of the system. It is easy to see that it is given by (pUjl .The (mass) density of the 
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star particles is similarly given by 

2Zo (A) Zo (B) + Zo (C) Zo (C) 



Since every star particle occupies five sites on the lattice, the number density of stars is 



(47) 



•'st.n 



0st/5. (48) 



The total partition function can be used to obtain the thermodynamics of the system 
with a calculation that is formally identical to the one carried on in the previous section for 
the case of monomeric particles and squares. The free energy per site is still given by (|26|). 
but Q2 is the polynomial 

Q2 = — + 2ah + c^. (49) 

As in the previous case, the calculations are done in the grand canonical ensemble and are 
carried out at constant chemical potentials. Hence, the above free energy represents l3PstVo 
where (3 is the inverse temperature, fn represents the volume of the unit cell of the lattice 
and Pst is the osmotic pressure across a membrane permeable to the star particles. 
The osmotic pressure Pg across a membrane permeable to the solvent particles is given by 
I3PsVq = (3PstVQ — In 77. If the solvents represent voids, then P^ represents the conventional 
pressure of the lattice system js^. 

In order to determine which phase is the stable one at some temperature, we must find 
the free energy of all the possible phases of the system as a function of w. 



2. /i > 

Let us first consider the case of positive chemical potential. As explained in the previous 
section, when fi is positive, the ground state of the lattice at zero temperature is represented 
by a pure solvent phase: all the lattice sites are occupied by monomeric species. As we start 
increasing the temperature the density of stars will increase. For positive values of /i, as in 
the previous case, only one phase is present and, therefore, no phase transition is observed. 
In this case, the solvent is percolating at every temperature and what we need to investigate 
is the percolation of stars. 

In order to study the percolation of the star particles, we introduce the probabilities 
Rm{c() < 1 as above. For example, Rm{C) < 1 is the probability that a site in a C state 
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at the mth generation is connected to a finite cluster of stars at higlier generations. At tlie 
1-cycIe fixed-point, each one of these Rm{(y) approaches its fixed-point value given by the 
solution of an equation of the form (|27p given earlier. By considering the configurations 
shown in Figures ^1 and ^1 and following the same steps shown in the previous section for 
the square particles, we can write: 

^ ' (a + w^s) ^ ' 

wh^ + (2w^i?(A) + w^)s''a + ?>w^R{kfsa^ + Ri.kfa' + /g(C)/g(B)^c6^ 
^ '~ w'^s^ + (2^2 + w^)s'^a + 3w2sa2 + + ch^ ' ^ ' 

[a + w^s) 

This system of equations does not depend explicitly on the chemical potential n but it 
depends on it implicitly. In fact, for every value of T, and hence of w, a different chemical 
potential corresponds to different values of s, a, h and c and consequently a different value 
for the i?(Q!)'s. 

Since R{q) represents the probability of the (occupied) site at the origin of the lattice 
being connected to a finite cluster of stars on one half of the lattice and since two half lattices 
are joined together at the origin to form the complete lattice, we define 

Pst = 1 - R{Cf (53) 

as the percolation probability for the stars. This quantity represents the probability to have 
the origin occupied by the core of a star connected to an infinite cluster of stars on both 
sides of the origin. As explained before 3^, there are many other possible definitions for 
the percolation probability. The one above is the one that seems the most proper choice for 
this particular case. The behavior of pst as a function of 0st for different (positive) values of 
/i, is shown in Figure |^ As expected, the percolation threshold increases as the chemical 
potential increases. The athermal case corresponds to /x = +oo,e such that fj,e is a 
positive constant. 

Figure |221 shows how the percolation threshold depends on the value of the chemical 
potential while Figure EH shows the dependence of the critical temperature. As expected, 
an increase in the chemical potential for the smaller particles makes it more and more difficult 
for the bigger particles to percolate. This explains the increase in the critical temperature 
and the percolation threshold as fj, increases. 
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For any given value of the chemical potential, the percolation threshold and the critical 
temperature are lower in the case of star particles as compared to the case of square parti- 
cles presented in the previous section. Consequently, our results describe the fact that, as 
observed in experimental results, an increase in the size disparity among the particles that 
form a physical system leads to a decrease in the percolation threshold. It is necessary to 
have a smaller fraction of sites occupied by the star particles, as compared to the square 
particles, in order to observe percolation. 

3. ii<0 

If negative values of fi are considered, the problem becomes more interesting and remark- 
ably more complicated, even compared to the problem of square and monomeric particles 
analyzed in the previous section. As explained in the previous section, a negative value for 
fi means that, according to thermodynamics, the system will be stable at low temperature 
if all the sites are covered by star particles. 

As it happens with the system formed by monomeric and square particles, what can be 
expected at low temperature, or at very low concentrations of solvent molecules, is an ordered 
structure of stars which disorders at high temperature, so that there is an abundance of star 
particles in the system. Thus, we could have a phase transition between an ordered phase, 
stable at low temperature and a disordered phase, stable at higher temperature. Figure ESI 
shows the free energy of the system in the case of fi = —1. 

By using the set of equations introduced above, the free energy that one obtains is the one 
represented by the open circles in the figure. This cannot be the free energy corresponding 
to the stable phase at every temperature for the same reason as given in the previous section, 
since it becomes negative below some critical temperature. 

In order to capture the ordered phase for this problem, it is necessary to have a calculation 
scheme that goes beyond the one introduced above. This is exactly what had to be done 
in the case of monomeric and square particles in the previous section. The main problem 
here is related to the fact that even though we had a somehow clear understanding of what 
the geometrical nature of the ground state should have been, its mathematical description 
is everything but simple. If it is necessary to solve one recursive equation in one variable 
there are very well known theorems that provide necessary and sufficient conditions for the 
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convergence of the recursion relation of one variable, and consequently for the solution of the 
problem. But, if we have a system of two or more equations in two or more variables, even 
as simple as possible but not linear, there are no conditions that guarantee the convergence 
of the recursion relations. Consequently, the solution is obtained through a process of trial 
and error. Unfortunately, even though a large number of different systems of equations have 
been solved, we have not been able to find a scheme that works for all the systems that 
we have considered. Even if we can immediately determine what approach will not work, 
for example in the case of the stars we cannot expect the one-cycle approach to be able to 
describe the ground state for the reasons stated above, we cannot tell a priori what recursion 
scheme will work. In the case of the system made of monomeric and star particles, a very 
complicated scheme has been used. The details of the equations used to solve the problem 
are not very enlightening. So the equations are not reported here. But here we will give the 
rationale behind the choice of this particular set of recursion relations for this problem. 

The key is to distinguish the directions along which one moves on the lattice and the 
position along which a star is placed inside a given plaquette. Figure |2E1 shows how this can 
be achieved. We distinguish the directions labeled 1, 2, 3 and 4 in the figure because of the 
very peculiar nature of the ground state. It is necessary to use this description in order to 
be able to capture the ground state that has a very regular pattern and where, as explained 
above in the case of the percolation of square particles, at low temperature when the lattice 
is almost completely covered by stars there are very strong correlations between sites that 
are second or third neighbors on the lattice. In this case, in order to get a solution it was 
necessary to introduce a symmetry breaking field to make the ground state unique. The 
symmetry breaking field makes some direction more favored than the others, thus creating 
a unique ground state for the system. There are four different fixed points that come out 
from the recursion relations. These different fixed points have to be carefully accounted for 
in order to obtain the proper free energy for this system. 

In this case, the approach is more complicated than in the case of the square particles. It 
is necessary to use a multi-cycle method. In this case, what is checked is the difference of the 
values of all the ratios at two levels that are n-generations apart on the lattice. Whenever 
<^m — Oim+n for all the ratios a is less than the desired tolerance we assume the fixed point of 
the solutions has been reached. Wherever the ordered phase exists (at low temperature) for 
the case of stars on a square Husimi lattice, we have a four-cycle structure of the fix-point. 
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so that the quantity am — ctm+4 is the one that is going to vanish. 

By using this second method, it is possible to obtain a second solution for the problem. 
The corresponding free energy curve is the one represented by the circles in Figure 031 The 
solutions obtained with the two different schemes coincide at high temperatures, but are 
different at low temperatures. The temperature at which the ordered phase appears and 
becomes the stable solution is, in the case of /i = —1, Tqd — 8.4. We use the subscript 
OD, as done before in the case of the square particles, to remind that this transition is a 
phase transition between an ordered and a disordered phase. It is possible to observe that 
the stable phase at every temperature has a positive free energy, as expected. 

The entropy corresponding to the free energy shown in Figure 123 is shown in Figure 
1771 As expected, below Tqd the entropy of the ordered phase is lower than the entropy of 
the metastable continuation of the disordered phase. The entropy of the metastable phase, 
though, drops much more rapidly, and goes to zero at a finite temperature. This temperature 
is the Kauzmann temperature of the system. 

The results for the percolation probability as a function of the temperature are shown 
in Figure EHl It is obvious that the strength of the percolation process is very different 
for the two phases. The formation of a percolating cluster of solvent molecules is much 
easier in the disordered phase than in the ordered one. There is a wide temperature range, 
between ~ 2.70 and ~ 7.45 in the case /i = —1, in which there is at least one percolating 
cluster in the disordered phase but no percolation occurs in the ordered phase. The region 
below 2.70, limited by the dashed line in Figure EHl does not represent physical states of 
the disordered phase because it is below the Kauzmann temperature and it corresponds 
to a negative entropy. This figure shows how the percolation threshold for the two phases 
is very different. In particular the solvent density must be at least ~ 0.32 in order to 
have percolation in the ordered phase while it is much lower in the metastable extension 
of the disordered phase. During the lowering of the temperature, if the phase transition is 
avoided, the system will percolate all the way down to the Kauzmann temperature. Figure 
EHl summarizes the results obtained as the value of the solvent chemical potential is changed. 
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V. CONCLUSIONS 



If we compare the results described above, we can draw some interesting conclusions 
about the effect of the size and shape of the particles on the percolation properties of these 
systems. 

The two systems that have been considered on the Husimi lattice have one common 
species, the monomeric one, and a different one, square particles in the first case and star 
particles in the second case. 

Both systems present common features. For instance, a metastable phase is present that 
is percolating very strongly in a temperature interval where the stable phase has a very 
weak percolating cluster (a much lower value of p) or even no percolating cluster at all. 
This observation suggests how it might be useful to prepare a system in a metastable state 
in some cases in which a percolating network of particles of some kind is needed in order 
to enhance one or more physical properties of the composite system. The presence of the 
star particles lowers the percolation threshold that decreases as the size disparity between 
different particles present in the system increases since the each star particle occupies more 
sites than a square particle. The introduction of the star particles, as compared to the 
square particles has an effect on the thermodynamics of the system as well. If we consider 
the case of negative values of the monomeric species chemical potential, we observe in both 
cases a continuous transition between an ordered phase that is stable at high temperature 
and a disordered phase that is instead stable above a critical temperature. The system that 
contains the star particles, though, has a much lower value of the critical temperature, as 
compared to the system containing the square particles. The presence of the larger star 
particles makes the disordering of the system easier as compared to the case in which the 
smaller square particles are present. So the system has a critical temperature that is lower 
in the case of the presence of star particles than in the case of the square particles. 

We have recently also looked at the percolation of stars and solvent particles on a Bethe 
lattice In this case, the most striking difference, as compared to the case of the Husimi 
lattice, is the different nature of the percolation transition since on the Bethe lattice we 
obtain a discontinuous transition for the percolation probability. At this stage of the re- 
search, the results obtained do not allow us to draw any conclusions about the reason of the 
different nature of the phase transition between the ordered and the disordered phase on the 
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two different lattice structures that we have used, the Husimi lattice and the Bethe lattice. 
We believe that the size disparity coupled with the nature of the lattice plays an important 
role in determining the nature of the transition. Further investigations are needed to give 
a complete answer to this question and they will be pursued in the future. In particular, it 
would be very interesting to study the behavior of larger and more complicated particles on 
both a Husimi lattice and a Bethe lattice. Of course, an increase in the complexity of the 
particles that are treated would likely lead to a tremendous increase in the complexity of 
the recursive equations that need to be solved. This would represent a problem mainly for 
the description of the ground state of the phase rich in larger particles. 

Now that we are able to describe the properties of systems made of particles of different 
sizes and shapes, we move forward and include the effect of having a polymer matrix in 
which the particles are embedded js^. 
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FIGURE CAPTIONS 
FIG. 1: Portion of an infinite lattice known as Husimi tree. 



FIG. 2: Possible states for the mth level site in the case of percolation of solvent and square 
particles. A, B and S refer to the three possible states of the base site of the plaquette, see text 
for details. 

FIG. 3: Percolation of squares and solvent molecules: possible configurations of the sites in the 
mth level square when the base is in the S state. 

FIG. 4: Percolation of squares and solvent molecules: possible configurations of the sites in the 
mth level square when the base is in the B state. 

FIG. 5: Percolation of squares and solvent molecules: possible configurations of the sites in the 
mth level square when the base is in the A state. 

FIG. 6: Probability that a square particle at the origin of the lattice is connected to an infinite 
cluster of square particles spanning the entire lattice as a function of the density of squares. 

FIG. 7: Probability that a square particle at the origin of the lattice is connected to an infinite 
cluster of square particles spanning the entire lattice as a function of the temperature of the system. 

FIG. 8: Percolation threshold as a function of the chemical potential for the percolation of squares 
on a Husimi tree filled with square and solvent particles. 

FIG. 9: Critical temperature as a function of the chemical potential for the percolation of squares 
on a Husimi tree filled with square and solvent particles. 

FIG. 10: Free energy as a function of temperature for = —1. See text for details. 
FIG. 11: Entropy as a function of temperature for /i = — 1. See text for details. 
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FIG. 12: Solvent percolation probability as a function of temperature for fi = —1 in the presence 
of square particles. The dotted line corresponds to the Kauzmann temperature. 

FIG. 13: Solvent percolation probability as a function of solvent density for ^ = —1 in the presence 
of square particles. 

FIG. 14: Critical temperature as a function of the chemical potential for the percolation of 
monomeric filler particles on a Husimi tree fihed with square and solvent particles. 

FIG. 15: Percolation threshold as a function of the chemical potential for the percolation of 
monomeric filler particles on a Husimi tree fihed with square and solvent particles. 

FIG. 16: Possible states for the mth level site in the case of percolation of solvent and star particles. 
A, B, C and S refer to the four possible states of the base site of the plaquette, see text for details.. 

FIG. 17: Percolation of stars and solvent molecules: possible configurations of the sites in the mth 
level square when the base is in the S state. 

FIG. 18: Percolation of stars and solvent molecules: possible configurations of the sites in the mth 
level square when the base is in the A state. 

FIG. 19: Percolation of stars and solvent molecules: possible configurations of the sites in the mth 
level square when the base is in the B state. 

FIG. 20: Percolation of stars and solvent molecules: possible configurations of the sites in the mth 
level square when the base is in the C state. 

FIG. 21: Probability that a star particle at the origin of the lattice is connected to an infinite 
cluster of star particles spanning the entire lattice as a function of the density of stars. 
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FIG. 22: Probability that a star particle at the origin of the lattice is connected to an infinite 
cluster of star particles spanning the entire lattice as a function of the temperature of the system. 

FIG. 23: Percolation threshold for star particles as a function of the chemical potential for the 
percolation of stars on a Husimi tree filled with star and solvent particles. 

FIG. 24: Critical temperature as a function of the chemical potential for the percolation of stars 
on a Husimi tree filled with star and solvent particles. 

FIG. 25: Monomeric and star particles: free energy as a function of temperature for = —1. 

FIG. 26: Labeling of the tree used to obtain the ground state of the system for the case of star 
and solvent particles on a Husimi tree with a ground state made of a pure star phase. 

FIG. 27: Monomeric and star particles: entropy as a function of temperature for fj, = —1. 

FIG. 28: Monomeric and star particles on a Husimi tree: solvent percolation probability as a 
function of temperature for /i = — 1. The dashed line corresponds to the Kauzmann temperature. 

FIG. 29: Monomeric and star particles on a Husimi tree: critical temperature as a function of the 
chemical potential for the percolation of monomeric filler particles on a Husimi tree. 



34 



m = 

FIG. 



35 




FIG. 3 



37 




FIG. 4 



38 




40 




T 



FIG. 7 



41 




FIG. 8 



42 




FIG. 9 



43 




44 




45 




46 



O one-cycle solution 
• two-cycle solution 



ffl 
O 
Ol 
O 

o 

o 
o 
o 
o 



ox 



0.2 0.3 



FIG. 13 



47 



disordered phase 
ordered phase 



-3 -2 



FIG. 14 



48 



• disordered phase 
■ ordered phase 



-4 -3 -2 
1^ 



-1 



FIG. 15 



49 



54 




FIG. 21 



55 




56 




FIG. 23 



57 




0.2 0.4 0.6 0.8 
1^ 



FIG. 24 



58 




FIG. 25 



59 




FIG. 26 



60 




61 




62 




FIG. 29 



63 



